Ground and Metastable States in 7-Ce from Correlated Band Theory 
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Multiple energy minima of the LDA+U energy functional are obtained for 7-Ce when it is 
implemented in a full potential, rotationally invariant scheme including spin-orbit coupling, and 
different starting local configurations are chosen. The lowest energy solution leads to a fully spin 
polarized 4f state and the lattice constant of 7-Ce. The higher energy local minima (additional 
self-consistent solutions) are shown to be strongly indicative of crystal electric field and multiplet 
excitations. 



The isostructural a - 7 phase transition in Ce, and the 
character of the phases separately, have been a subject 
of continious theoretical efforts for decades. Experimen- 
tally, the low temperature fee a phase of Ce transforms 
to the fee 7 phase at the high temperatures with a large 
increase of the volume. When pressure is applied, the 
crystal collapses back into the a phase in a first-order 
phase transition. The room temperature 7 — a phase 
transition occurs at ~ 7 kbar pressure with the decrease 
of the volume of « 15%. 

Several theoretical models have been suggested and 
they are described in some detail in Rcf. (lj. The differ- 
ence between the models arises from the different treat- 
ments for the Ce 4f electrons. The promotional model || 
assumed depopulation of the 4f level upon compression, 
but does not agree with the results of various experiments 
which indicate little change in the 4f occupation. Band 
structure calculations also indicate that promotion is too 
high in energy to drive the transition, and the 4f occupa- 
tion remains near unity. || The Mott transition model 
H assumes the transformation of localized to band-like 
4f state with the decrease of volume and is similar to 
the Mott-Hubbard metal-insulator transition. First prin- 
ciples calculations using the self interaction correction 
(SIC) to the LDA || produce total energy minima for 
localized (7) and band (a) states with formally the same 
SIC total energy variational functional, and are in ac- 
cord with the Mott transition model. However, the use of 
atomic sphere approximation (ASA) , which sphericalizes 
the potential and density within large atomic spheres, 
limits the ability of these calculations to produce quanti- 
tative total energy description of the system with highly 
non-spherically symmetric 4f electron charge/spin densi- 
ties. 

The other viable theoretical model is the Kondo vol- 
ume collapse (KVC) model ](J. The essence of the KVC 
is an assumption of localized f states in the both a and 7 
phases of Ce. The a phase consists of a mixed valence 4f 
state while the 7 phase has almost integer 4f occupancy. 
The phase transition is due to the entropy contribution 
of the localized 4f state. The relation between KVC and 
Mott transition model was analyzed recently in Ref. |l]] 
and quantitative arguments in favor of KVC model were 
provided. 



The local (spin) density approximation (LDA) has had 
tremendous success in the quantitative description of a 
wide varietly of solids, (7j but the rare earths are ex- 
treme cases where the LDA description is inadequate. 
The "correlated band theory" (LDA+U) approach has 
had much success in the treatment of correlated mag- 
netic insulators where LDA results are incorrect. In this 
paper we explore not the 7 — a transition per se, but 
the 7 phase itself, by an application of the LDA+U ap- 
proach that includes the new features (1) a full potential 
method is used, which enables the novel features we un- 
cover, and (2) the LDA+U method is applied to a metal 
where the interaction of the local orbitals with the con- 
duction states is a central part of the physics. For the 
first time, multiple LDA+U energy minima (ground plus 
metastable states) are obtained within a specified bro- 
ken symmetry, and we discuss how the various states are 
related to 4f excitations in the 7 phase. 

The LDA+U method || is based on an energy which 
is a functional of the spin densities {p s , s = ±1}, and the 
occupation matrices {n s m m ,} of the 4f orbitals labelled 
by their azimuthal projections m = mi. There is every 
reason to expect that, within the allowed values of p s and 
n s , there can be local energy minima as well as the abso- 
lute minimum, which represents the ground state energy 
of the system. (Such local minima are rarely found within 
LDA except in magnetic systems.) The formal meaning 
of these local minima, as well as the formal underpinnings 
of the LDA+U approach itself, remain to be settled, but 
- like the Kohn-Sham eigenvalues (the band structure) 
which have little formal meaning but immense practical 
importance - these minima will be shown to bear a close 
relationship to local excitations of the 4f shell. 

Our results are based upon the full-potential linearized 
augmented plane wave method (LAPW) as the ba- 
sis for total energy calculations with the rotationally in- 
variant LDA+U functional fiofl , with spin-orbit coupling 
(SOC) is included self-consistently |jyj. Literature values 
fL2f are used for the on-site repulsion U = 6.1 eV and 
exchange J = 0.7 eV (Slater integrals, F$ = 6.10 eV, 
F 2 = 8.34 eV, F 4 = 5.57 eV, F 6 = 4.12 eV ). Ferromag- 
netic spin alignment is assumed (allowed, not imposed) 
since local moments are present in paramagnetic high 
temperature 7 phase of Ce [n3| . 
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The key to finding the various solutions (local minima) 
is to start from atomic densities and different 'guesses' 
for the 4f occupation matrices n s and then obtain self- 
consistently both spin/charge densities and occupation 
matrices. As an intuitive guess for 4f occupations we use 
the fully spin-polarized state (Tr n'=l, Tr n^=0) and 
choose various orbital characters for nJ m m , . 

Ground State. The calculated equilibrium lattice con- 
stant (cf. Table I) is very close to the experimental value 
of 7-Ce at zero pressure and room temperature fl4f| . The 
calculated bulk modulus is about 25% larger than the ex- 
perimental value, similar to other calculations that have 
treated the localized character of the 4f state in 7-Ce. 

For the 7 phase volume (lattice constant a=9.76 a.u.) 
we find the lowest solution to be fully spin polarized with 
the f occupation of 1.04. This state is primarily m;=-2 
(rrij — — |), with some mixing in of the spin-majority 
mi = 2 state. The LDA+U yields significant enhance- 
ment for the absolute values of both spin and orbital 
magnetic moments (anti-aligned in accord with the 3rd 
Hund's rule) in comparison with LDA results. 

The electronic density of states for the 7-Ce ground 
state is shown in Fig. 1 in comparison with the re- 
sult of relativistic (with SOC) LDA calculations. The 
f-majority peak at the bottom of valence band (ss 2.5 eV 
below the Fermi level) indicates the position of the lo- 
calized 4f state, in quantitative agreement with the res- 
onant photoemission measurements Jl5[ . The difference 
between LDA+U and LDA is due to the localization of 
the 4f state, which removes both majority and minority 
4f states from the vicinity of the Fermi level. Since the 
4f state is not well separated from the valence band (cf. 
Fig. 1) we conclude that in spite of its localized character 
the 4f state in 7-Ce cannot be treated as core-like. 

Metastable States. As mentioned above, several self- 
consistent solutions corresponding to different (local) 
minima of the LDA+U energy functional are possible. 
Different minima must be searched for by exploring var- 
ious regions of the underlying space. Our initial starting 
points consisted of fully spin polarized 4f x with differ- 
ent orbital characters {mi}. Our initial and final (self- 
consistent) majority spin occupation matrix elements are 
shown in Table II. In the LDA+U calculations the spin 
polarization enters as an effective magnetic field provid- 
ing the spin quantization axis and the direction of the 
quantization axis is chosen along z([001]). 

The LDA+U minima can be considered as mean-ficld- 
like projections (in the Fock space) of many-body wave- 
functions on the single-particle angular basis set. The 
interpretation is much simpler for an f 1 ion than it would 
be for a multi-electron ion because the {mi,m s } and 
{ J, J z } representations are unitarily related. In the sim- 
ple crystal-electric-field (CEF) model for the paramag- 
netic Ce 3+ ion (l7) the ground state is formed by one of 
the T7 and Tg' 2 doublets from J = 5/2 multiplet In 



order to compare the CEF model |17| with the LDA+U 
results, we apply to the LDA+U solutions the unitary 
transformation from {mi,m s } representation to {J, J z } 
representation with J = 5/2, 7/2. The states from Table 
III are then classified as follows: 

Solution (1) - Ground State: It has 96% mi — —2, with 
3% mi = 2 mixed in. The transformation to {J, J z } 
representation yields 69% of 1 5/2, — 3/2 > and 1% of 
|5/2, 5/2 > states from J = 5/2 multiplet, plus 27% of 
|7/2, -3/2 > and 3% of |7/2, 5/2 > states from J = 7/2. 
We conclude that the Solution 1 consists of 70% states 
from the J = 5/2 multiplet which are the linear combi- 
nation of the CEF levels T7 and Y\, and 30% states from 
the J — 7/2 multiplet. 

Solution (2): Roughly equal amounts of m; = —3 and 
mi = +1) in this solution indicates coupled mj = —5/2 
and mj — 3/2 states. It consists of 55% states from 
J = 5/2 multiplet (combinations of the CEF levels Tj 
and Tg) and 45% of states from J = 7/2. 
Solution (3), involving admixture of m; = — 1 with 25% 
mi = +3, which transforms to an admixture of 42% of 
|5/2, -1/2 > >) state with the 58% of |7/2, 7/2 >. 
Solution (4), pure mi = 0, corresponds to the combina- 
tion of 42% of |5/2, 1/2 > (|r| >) state with the 58% of 
|7/2, 1/2 >, and lies 232 meV above the ground state. 

These correlated band structure calculations are pro- 
ducing metastable solutions that contain the CEF states 
from the conventionally presumed lowest multiplet for 
Ce 3+ . In addition, there is a considerable fraction of the 
states from the first excited J7/2 multiplet that must be 
understood. 

There are two important reasons why our solutions do 
not correspond directly to CEF levels as normally con- 
sidered. First, the CEF picture assumes ideal cubic sym- 
metry, whereas a 4f state in 7-Ce is actually surrounded 
by twelve atoms whose own moments are oriented ran- 
domly, hence breaking cubic symmetry. In addition the 
non-cubic nature of the 4f density perturbes the conduc- 
tion electron density, which leads to a non-cubic local 
field. Secondly, in our calculations we have artificially or- 
dered the spins of the magnetic ions and allowed orbital 
moments, which reduces out site symmetry to tetragonal. 

The comparison between LDA+U calculated splitting 
scheme with the results of inelastic neutron scattering 
experiments fl9| is shown in Fig. 2. It is clearly seen 
that Solutions 1-3 lie well within the range of low-energy 
excitations peak and reflect the mixed CEF and spin- 
orbit excitations. The energy position of Solution 4 (232 
meV) agrees quantitatively with the experimental inelas- 
tic peak at 260 meV. This peak is usually interpreted 
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calculations suggest that this interpretation is oversim- 
plified since the states from both lower 2 F 5 / 2 and first 
excited 2 FV/2 multiplets are mixed in the ground and ex- 
cited states of 7-Ce. The spectrum from present LDA+U 
calculations are in better agreement with the experiment 
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[19 1 than previously reported results of SIC calculations 
[2Cfl (86-100 meV and 130 meV). 

In order to classify the LDA+U solutions we transform 
them from the compex to the cubic spherical harmonics. 
Solution (1) involves admixture of xyz and (x 2 — y 2 )z 
states; Solutions (2) and (3) consist of mixed {x 2 — z 2 )y 
and (y 2 — z 2 )x states; and Solution (4) has (5z 2 — 3r 2 )z 
symmetry. All the LDA+U solutions have the tetrago- 
nal symmetry assumed in the calculations. The local- 
ized states without SOC are formed by CEF states in 
{L, L z } representation and for L = 3 these are L2, L4, T5 
states Iplfl . In fully spin-polarized case the spin-orbit 
coupling £(/ • s) is reduced to £,(l z s z ) and the eigenval- 
ues of localized "CEF+SOC" model Hamiltonian are eas- 
ily obtained. The "CEF+SOC" coupling scheme is pre- 
sented in Fig. 3 and shows how the LDA+U solutions 
result from CEF eigenstates coupled by spin-orbit cou- 
pling. From the total energy differences between LDA+U 
solutions we then derive the CEF splitting parameters 
A 2:5 = T 2 - T 5 = 90 meV and A 4 , 2 = T 4 - T 2 = 107 
meV and SOC parameter £ of 66 meV (SOC splitting = 
(7/2)£ = 231 meV). 

The scheme in Fig. 3 allows more eigenstates than the 
calculated local minima of the LDA+U total energy func- 
tional. In each case the occupation matrices n mi . m j for 
these additional states have the non-zero elements in the 
same {mi, m^} sub-space as one of the metastable states. 
As a result the minimization procedure yields the lowest 
total energy state in given {mi,m[} sub-space and the 
higher states do not appear among variational LDA+U 
solutions. 

The physical origin of the CEF splitting can be under- 
stood to be due to the anisotropy of the mixing interac- 
tion between conduction band and localized f-states [^2| . 
We obtain the CEF splitting for 7-Ce of the order of 100 
meV comparable with the SOC splitting of 231 meV. In 
both ground and excited states the CEF levels are then 
mixed by SOC yielding low-energy excitations of the or- 
der of 20-50 meV and high-energy excitations about 230 
meV. 

The results of LDA+U calculations allow us to con- 
clude that the interaction between conduction band and 
localized f-states plays an important role for the both 
ground and excited states in 7— Ce. As a result, the Ce 
f-states are not core-like and should be treated in terms of 
quantum impurity (Anderson) model |2^| rather than as 
localized states of Ce 3+ ion in the cubic CEF. The present 
LDA+U calculations correspond to the numerical solu- 
tion of the lattice Anderson model in static mean-field 
approximation with the assumed ferromagnetic order. 
This simplification does not allow us to describe quantita- 
tively the paramagnetic (disordered local moment) high- 
temperature state of 7— Ce observed experimentally. The 
treatment of this paramagnetic state requires the use of 
dynamic mean-field parametrization for the self-energy 



Recently, Solovyev et al. [|23j proposed a correction to 
the LDA+U total energy functional [|ij when SOC is 
included and non-collinear magnetic configurations are 
considered. This correction accounts for additional con- 
tributions to the exchange energy due to non-zero spin- 
off-diagonal elements of the occupation matrix {n^m'}- 
We performed LDA+U calculations with this spin-off- 
diagonal correction |23) and found that the spin-off- 
diagonal occupations are very small and have minute ef- 
fects on the values of spin and orbital magnetic moments, 
and the total energies (less than 2 meV) for ground and 
metastable states (cf. Table II) in 7— Ce. 

To summarize, we have obtained the ground and three 
metastable states from relativistic (with spin-orbit cou- 
pling) spin-polarized full-potential LDA+U calculations 
for fee Ce. The ground state has equilibrium lattice con- 
stant of 7-Ce. Our calculations reproduce the localized 
character of 4f states in 7-Ce, which however cannot be 
treated as a part of atomic core. Analysis of various 
LDA+U solutions allows us to make a comparison be- 
tween correlated band theory results and CEF model. 
The excitation energies calculated from the total energy 
differences between LDA+U solutions are in reasonable 
agreement with the experimental data. 
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FIG. 1. Spin-resolved total and partial 4f densities of states 
for ordered 7-Ce. Majority are plotted upwards, minority 
downwards. The 4f partial DOS from LDA is dotted. The full 
line gives the LDA+U DOS, with the 4f contribution filled in. 
The 4f DOS peaks near the Fermi level in both spins are well 
removed from the Fermi level by LDA+U. 
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FIG. 2. Schematic sketch of the experimental inelastic neu- 
tron scattering intensity in comparison with the energy posi- 
tions of LDA+U calculated ground and excited states [l]-[4]. 
Additional inferred localized excitations are shown (dotted 
lines) . 



CEF + SOC 




L3J = (x 2 - z 2 )y + (y 2 - 2? )x 

L2J = (x 2 - z-)y +(y-- z 2 )x 
[l] = xyz + (x 2 - y 2 )z 



FIG. 3. A schematic sketch for the relation of the CEF 
levels to the LDA+U calculated ground and excited states 
[l]-[4] (full lines). Additional possible localized excitations 
are shown (dotted lines). 



TABLE I. The equilibrium lattice constant (a, a.u.) and 
bulk moduli (B, kbar) as a result of LDA+U calculations. 



LAPW 


LMTO-ASA | 


} LMTO-ASA (T| 


Exp. § 


LDA+U 


SIC 


f-core LDA f-core GGA 




a 9.83 


9.58 


9.69 10.02 


9.76 


B 296 


310 


312 288 


210, 244 



TABLE II. The elements of the occupation matrix n m , m / 
(spin - majority) (initial assignment n° and self-consistent 
n sc f)\ spin (2 < S z >), orbital (< L z >) moments and 
the total energy increase with respect to the ground state 
(AE, meV) for different LDA+U self-consistent solutions in 
the order of increasing total energy. The only those ele- 
ments of occupation matrix which are bigger than 0.01 are 
shown. The spin-minority occupation matrix is almost zero 
since there is a complete spin-polarization of the 4f-shell. 



n° m 


m' 


sc f 

n , 




2 < S z > 


< L z > 


AE 


Solution 1 


m, m! 


-2 


2 m,m' -2 


2 


1.18 


-1.87 





-2 


1 


-2 0.960 


0.143 








2 





2 0.143 


0.030 








Solution 2 


m, m! 


-3 


1 m,m' -3 


1 


1.18 


-0.75 


19 


-3 


1 


-3 0.438 


-0.484 








1 





1 -0.484 


0.552 








Solution 3 


m, m! 


-1 


3 m,m' -1 


3 


1.18 


0.001 


52 


-1 


1 


-1 0.737 


-0.423 








3 





3 -0.423 


0.253 








Solution 4 


m, m! 





m, m! 




1.14 


-0.004 


232 





1 


0.979 
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